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Abstract Using the 32-m Medicina, 45-m Nobeyama, and 100-m Effelsberg tele- 
scopes we found a statistically significant velocity offset AV w 27± 3 m s^' (Iff) 
between the inversion transition in NH3(1,1) and low -J rotational transitions in 
N2H+(l-0) and HC3N(2-1) arising in cold and dense molecular cores in the Milky 
Way. Systematic shifts of the line centers caused by turbulent motions and velocity 
gradients, possible non-thermal hyperfine structure populations, pressure and optical 
depth effects are shown to be lower than or about 1 m s^' and thus can be neglected 
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in the total error budget. The reproducibihty of AV at the same facihty (Effelsberg 
telescope) on a year-to-year basis is found to be very good. Since the frequencies 
of the inversion and rotational transitions have different sensitivities to variations in 
jj. = me/nip, the revealed non-zero 4 V may imply that jj. changes when measured at 
high (terrestrial) and low (interstellar) matter densities as predicted by chameleon- 
like scalar field models - candidates to the dark energy carrier. Thus we are testing 
whether scalar field models have chameleon-type interactions with ordinary matter 
The measured velocity offset corresponds to the ratio Ap,/ii = (|lspace — Miab)/ Miab 
of (26 ±3) X 10-*^ (la). 



1 Introduction 

This contribution sums up our results of differential measurements of the electron- 
to-proton mass ratio, pL = me/nip, carried out at the 32-m Medicina, 45-m Nobeyama, 
and 100-m Effelsberg telescopes ll24l l29l l25l l26l . With high spectral resolution 
(FWHM ^30-40 m s^' ), we observed narrow emission lines (FWHM < 200 m s^' ) 
of N-bearing molecules arising in cold and dense molecular cores of starless molec- 
ular clouds located in the Milky Way. The cores are characterized by low kinetic 
temperatures T^m ^ 10 K, gas densities n ^ 10^ — 10^ cm^^ , magnetic fields Z? < 10 
jxG, and ionization degrees Xe ^ 10^^ H). The objective of this study is to probe 
the value of Aji/ji = (/ispace ^ Miab)/Miab which is predicted to be variable when 
measured at high (laboratory) and low (space) matter density environments lf33]| . 

A hypothetical variability of jj. is thought to be due to the scalar fields - candi- 
dates to the dark energy carrier, - which are ultra-light in cosmic vacuum but possess 
an effectively large mass locally when they are coupled to ordinary matter by the so- 
called chameleon mechanism |fT6| . Several possibilities to detect chameleons were 
discussed in ||2]|3. First laboratory experiments constraining these models have 
been recently carried out at Fermilab B2l and in the Lawrence Livermore National 
Laboratory 1381 . 

A subclass of chameleon models considers the couplings of a scalar field to mat- 
ter much stronger than gravitational and predicts that fundamental physical quanti- 
ties such as elementary particle masses may depend on the local matter density, p, 
||33]| . This means that a nonzero value of Ajx/jX is to be expected for all interstellar 
clouds irrespective of their position and local matter density because the difference 
Ap between the terrestrial environment in laboratory measurements and dense in- 
terstellar molecular clouds is extremely large, PQ/pcioud > lO"^. 

In the standard model (SM) of particle physics the dimensionless mass ratio fi = 
m^/mp defines the ratio of the electroweak scale to the strong scale since the mass of 
the electron is proportional to the Higgs vacuum expectation value and the mass 
of the proton m-g is proportional to the quantum chromodynamics scale Aqcd 13 . 
The SM is extremely successful in explaining laboratory physics, but it has serious 
problems in astrophysics where it completely fails to explain dark matter and dark 
energy. There are many extensions of the SM including supersymmetry and different 
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multidimensional theories which introduce new particles as possible candidates for 
the dark matter and additional scalar fields to describe the nature of dark energy. 

A concept of dark energy with negative pressure appeared in physics long be- 
fore the discovery of the accelerating universe through observations of nearby and 
distant Type la Supernovae 13611371 . Early examples of dark energy in the form of 
a scalar field with a self-interaction potential can be found in reviews [Sl and 1351 . 
If masses of the elementary particles are affected by scalar fields, one can probe 
the dimensionless constant ji at different physical conditions by means of high pre- 
cision spectral observations. Namely, since the inversion and rotational molecular 
transitions have different sensitivities to variations in nonzero Ajx/jj. causes 

an offset between the radial velocities {AV = Vmt — Vim) of co-spatially distributed 
molecules, which, in turn, provide a measure of zi/i//i. When the inversion line 
belongs to NHswe will herein call this procedure the ammonia method. 



2 The ammonia method 

NH3 is a molecule whose inversion frequencies are very sensitive to changes of /i 
because of the quantum mechanical tunneling of the N atom through the plane of the 
H atoms. The inversion vibrational mode of NH3 is described by a double-well po- 
tential, the first two vibrational levels lying below the barrier. The quantum mechan- 
ical tunneling splits these two levels into inversion doublets providing a transition 
frequency that falls into the microwave range |[T4l . 

The sensitivity coefficient to -variation of the NH3 {J,K) = (1,1) inversion 
transition at 24 GHz is ginv ~ 4.46 l9|. This means that the inversion frequency 
scales as 

(zito/co)inv = (<a - (o)/(o = AA6{Ajx/pL) , (1) 

where (O and d) are the frequencies corresponding to the laboratory value of and 
to an altered in a low-density environment, respectively. In other words, the inver- 
sion transition sensitivity to /i -variation is 4.46 times higher than that of molecular 
rotational transitions, where grot = 1 and thus, 

{A(0/(0)rol^All/ll. (2) 

In astronomical spectra, any frequency shift zi o is related to the radial velocity 
shift AVr (Yr is the line-of-sight projection of the velocity vector) 

AVr/c = [Vr - Vb)/c = (OJlab - COobs)/fiJlab , (3) 

where Vq is the reference radial velocity, c is the speed of light and Oiab, Oobs are 
the laboratory and observed frequencies, respectively. Therefore, by comparing the 
apparent radial velocity Vinv for the NH3 inversion transition with the apparent radial 
velocity Vrot for rotational lines originating in the same molecular cloud and moving 
with radial velocity Vo with respect to the local standard of rest, one can find from 
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Eqs. dll-S 

A 11/ 11^ 0.289(Vrot- Vinv)/c = 0.2894 V/c . (4) 
The velocity offset AV in Eq.© can be expressed as the sum of two components 

AV = AV^+AV„, (5) 

where AV^ is the shift due to /i -variation and AV„ is a random component caused 
by the inhomogeneous distribution of molecules, turbulent motions and velocity 
gradients, possible non-thermal hyperfine structure populations, pressure and optical 
depth effects, and instrumental imperfections (the so-called Doppler noise). 

We will assume that the Doppler noise component has a zero mean, (AVn) = 
km s^' , and a finite variance, Var{AV„) < °°. Then, the signal AVfj can be estimated 
statistically by averaging over a data sample 

(AV) = (AV^) , (6) 

and 

Var{AV) = Var{AV^) + Var{AV„) . (7) 

Equation (|7]i shows that the measurability of the signal (4V^) depends critically on 
the value of the Doppler noise, Var{AVn). The Doppler noise can be reduced by the 
appropriate choice of molecular lines and molecular clouds. 

To have similar Doppler velocity shifts the chosen molecular transitions should 
share as far as possible the same volume elements. NH3 inversion transitions are 
usually detected in dense molecular cores (« ~ 10^ — 10^ cm^^). Mapping of these 
cores in different molecular lines shows that there is a good correlation between 
the NH3, N2H+, and HC3N distributions HO] [El ET] [M) . However, in some clouds 
NH3 is not traced by HC3N. The most striking case is the dark cloud TMC-1, where 
peaks of line emissions are offset by 7 arcmin lf32l . Such a chemical differentiation 
together with velocity gradients within molecular cores is the main source of the 
unavoidable Doppler noise in Eq.©. Additional scatter in the AV„ values may be 
caused by the different optical depths of the hyperfine structure transitions and by 
external electric and magnetic fields discussed below. 



3 Selection of targets 

To minimize the Doppler noise component in Eq.(|5]l, a preliminary selection of 
molecular cores suitable for the most precise measurements of the velocity offsets 
AV between rotational and inversion transitions is required. In general, molecular 
cores are not ideal spheres and when observed at high angular resolutions they fre- 
quently exhibit complex substructures. The line profiles may be asymmetric because 
of non-thermal bulk motions. Taking this into account, we formulate the following 
selection criteria: 
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1 . The line profiles are symmetric and well described by a single-component Gaus- 
sian model. This selection increases the accuracy of the line center measurement. 
Multiple line components may shift the barycenter and affect the velocity differ- 
ence between molecular transitions because, e.g., the ratio NH3/HC3N may vary 
from one component to another. 

2. The line widths do not exceed greatly the Doppler width caused by the ther- 
mal motion of material, i.e., the non-thermal components (infall, outflow, tidal 
flow, turbulence) do not dominate the line broadening. This ensures that selected 
molecular lines correspond to the same kinetic temperature and arise co-spatially. 
For the molecules in question we require the ratio of the Doppler ^-parameters 
to be j3 = fe(NH3)/fe(HC3N) or j3 = fe(NH3)/fe(N2H+), be j3 > 1. 

3. The spectral lines are sufficiently narrow (b ^O.l~ 0.2 km s^ ' ) for the hyperfine 
structure (hfs) components to be resolved. This allows us to validate the measured 
radial velocity by means of different hfs lines of the same molecular transition. 

4. The spectral lines are not heavily saturated and their profiles are not affected by 
optical depth effects. The total optical depth of the NH3 hf transitions is T < 10, 
i.e., the optical depth of the strongest hf component is < 1. 

The kinetic temperature Tkin, and the nonthermal (turbulent) velocity dispersion, 
<Ttuib> can be estimated from the line broadening Doppler parameters b = V2<7v, 
of, e.g., the NH3 (1,1) and HC3N (2-1) lines. Here d,, is the line of sight velocity 
dispersion of the molecular gas within a given cloud. If the two molecular transitions 
trace the same material and have the same nonthermal velocity component, (7,, is the 
quadrature sum of the thermal Cth and turbulent Oiurb velocity dispersions. In this 
case a lighter molecule with a mass m/ should have a line width wider than with a 
heavier molecule with m/, > nii. 

For thermally dominated line widths (j3 > 1) and co-spatially distributed species 
we obtain the following relations (e.g., [[TOl ): 

Tidn = mmh{(yf - (yl)/k{mh - mi) , (8) 

and 

^''turb = i"^h (^h - )/{m-mi), (9) 

where k is Boltzmann's constant, and the thermal velocity dispersion C7th is given by 

(Tth,, = (fcrkin/m,)'/^ (10) 

It should be noted that HC3N is usually distributed in a volume of the molecular core 
larger than NH3 since N-bearing molecules trace the inner core, whereas C-bearing 
molecules occupy the outer part fSl. Such a differentiation may cause a larger non- 
thermal component in the velocity distribution of HC3N. If both molecules are 
shielded from external incident radiation, and the gas temperature mainly comes 
from heating by cosmic rays, then a formal application of Eqs.® and ^ to the 
apparent line widths provides a lower limit on Ttin and an upper limit on Cturb- In 
molecular cores where the only source of heating are the cosmic rays and the cooling 
comes from the line radiation mainly from CO, a lower bound on the kinetic temper- 
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ature is about 8 K IfTTII . Thus, point 2 of the selection criteria requires ^ 8-10 K 
if both lighter and heavier molecules are distributed co-spatially. 



4 Preliminary results 

In our preliminary single-pointing observations flSl we studied 41 molecular cores 
along 55 Hnes of sight. The NH3(1,1), HC3N(2-1), and N2H+(l-0) transitions were 
observed with the 100-m Effelsberg, 32-m Medicina, and 45-m Nobeyama tele- 
scopes. The analysis of the total sample revealed large systematic shifts and 'heavy 
tails' of the AV probability distribution function resulting in a poor concordance be- 
tween three mean values: the weighted mean (4y)w = 27.4±4.4 m s^' (weights 
inversionally proportional to the variances), the robust redescending M-estimate 
{AV)m = 14.1 ±4.0 m s^' , and the median ZiVmed = 17 m s^' . The scatter in the 
AV values reflects effects related to the gas kinematics and the chemical segregation 
of one molecule with respect to the other (approximately 50% of the targets showed 
J3<1). 

If we now consider molecular lines that fulfill the selection criteria, the sam- 
ple size is n = 23, i.e. two times smaller than the total data set. Nevertheless, for 
this reduced sample we obtained better concordance: (4y)w = 20.7 ± 3.0 m s^' , 
{AV)m = 21.5 ±2.8 m s"' , and AV,-ned = 22 m s"' . 

In this preliminary study we found 7 sources with thermally dominated motions 
which provide (ziy)w = 21.1 ± 1.3 m s"^ (ziy)M = 21.2 ± 1.8 m s"^ and ZiVmed = 
22 m s^' . 

The most accurate result was obtained with the Effelsberg 100-m telescope: 
AV = 27±4stat±3sys m s^' (the value is slightly corrected in |26|). The systematic 
error in this estimate is due to uncertainties of the rest frequencies of the HC3N(2- 
1) hfs transitions since uncertainties of the NH3(1,1) hfs transitions are less than 1 
m s^ ' 1211 . When interpreted in terms of the electron-to-proton mass ratio variation, 
this gives A^/^ = (26±4stat±3sys) x 10"^. 

Thus, for the first time we obtained an astronomical spectroscopic estimate of the 
relative change in the fundamental physical constant at the level of 10^^ which is 
10^ times more sensitive than the upper limits on Ajj./ jj. obtained by the ammonia 
method from extragalactic observations 191 l30l l28l [T3l . 



5 Mapping of cold molecular cores in NH3 and HC3N lines 

In this section we report on new observations in which we measure AV at different 
positions across individual clouds in order to test the reproducibility of the veloc- 
ity offsets in the presence of large-scale velocity gradients l26ll . From the list of 
molecular cores observed with the Effelsberg 100-m telescope, we selected sev- 
eral objects with symmetric profiles of the NH3(1,1) and HC3N(2-1) hfs transi- 
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tions. In these objects the Hne widths are thermally dominated, i.e., the parameter 
j3 = (7,,(NH3)/av(HC3N) > 1. The chosen targets are the molecular cores L1498, 
L1512, L1517B, and L1400K, which have already been extensively studied in many 
molecular lines H] ED |23l gl] gO] |6l . 

The inversion line of ammonia NH3(1,1) at 23.694 GHz and the rotation line 
of cyanoacetylene HC3N(2-1) at 18.196 GHz were measured with a K-band HEMT 
(high electron mobility transistor) dual channel receiver, yielding spectra with an an- 
gular resolution of HPBW ^ 40" in two orthogonally oriented linear polarizations. 
The measurements were carried out in frequency-switching mode using a frequency 
throw of 5 MHz. The backend was an FFTS (Fast Fourier Transform Spectrometer) 
operated with its minimum bandwidth of 20 MHz providing simultaneously 16 384 
channels for each polarization. The resulting channel separations are 15.4 m s^' for 
NH3 and 20. 1 m s^ ' for HC3N. We note, however, that the true velocity resolution is 
about two times lower, FWHM 30 m s^ ' and 40 m s^ ' , respectively ifTTl . The sky 
frequencies were reset at the onset of each scan and the Doppler tracking was used 
continuously to track Doppler shifts during the observations. Observations started 
by measuring the continuum emission of calibration sources and continued by per- 
forming pointing measurements toward a source close to the spectroscopic target. 
The calibration is estimated to be accurate to ±15% and the pointing accuracy to be 
higher than 10". 

After corrections for the rounded frequencies, the individual exposures were 
co-added to increase the signal-to-noise ratio S/N. The spectra were folded to re- 
move the effects of the frequency switch, and base lines were determined for each 
spectrum. The resolved hfs components show no kinematic substructure and con- 
sist of an apparently symmetric peak profile without broadened line wings or self- 
absorption features. The line parameters, such as the total optical depth in the tran- 
sition T (i.e., the peak optical depth if all hyperfine components were placed at the 
same velocity), the radial velocity V/sr, the line broadening Doppler parameter b, 
and the amplitude A, were obtained through fitting of the one-component Gaussian 
model to the observed spectra: 

r(v)=A.[l-expH(v))], (11) 

with 

k 

t{v) = T • ^ fl, exp [-(v - V, - , (12) 

1=1 

which transforms for optically thin transitions into 

k 

r(y)=A'.^fl,-exp[-(v-y,--y,,,-)V^'] • (13) 
1=1 

The sum in (fT2l l and ( fT3] ) runs over the A: = 18 and k = 6 hfs components of the 
NH3(1,1) and HC3N(2-1) transitions, respectively. 

To test possible optical depth effects, we calculated two sets of the fitting param- 
eters: (1) based on the analysis of only optically thin satellite lines with AF\ ^ 0, 
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and (2) obtained from the fit to the entirety of the NH3(1,1) spectrum including the 
main transitions with AFi =0, which have optical depths T « 1-2, as can be inferred 
from the relative intensities of the hfs components. The resulting Vis,- values for all 
targets were within the la uncertainty intervals. 

To control another source of errors caused by instrumental imperfections, we 
carried out repeated observations at 10 offset positions: two in L1498, L1517B, and 
L1400K, respectively, and four in L1512. The AV dispersion resulting from these 
repeated measurements is a{AV) = 2.0 m s^' . This dispersion is lower than the 
calculated 1(7 errors of the individual AV values, which ensures that we are not 
missing any significant instrumental errors at the level of a few m s^' . 

We also checked the velocity offsets AV obtained in our observations with 
the 100-m Effelsberg telescope in Feb 2009 and Jan 2010. The reproducibility of 
the velocity offsets at the same facility on a year-to-year base is very good for 
L1498, L1512 and L1517B (concordance within 1(7 uncertainty intervals), except 
for L1400K where the central point is probably an outlier. 

The comparison of the velocity dispersions determined from the NH3(1,1) and 
HC3N(2-1) lines did not show any significant variations with position within each 
molecular core. All data are consistent with thermally dominated line broadening. 
In particular, the following weighted mean values were determined: )3li498 = 1.24 ± 
0.03, j3Li5i2 = 1.32±0.05, jSusns = 1.11 ±0.01, and j3u4ooK = 1.23±0.04. The 
weighted mean values of the velocity dispersions for NH3 range between C7li512 = 
78 ± 1 m s^' and (7li4ook = 86 ± 1 m s^' , and for HC3N between Ousii = 59 ± 2 
m s~' and Ousm = 74 ± 1 m s^' . This can be compared with the speed of sound 
inside a thermally dominated region of a cold molecular core that is defined as (e.g., 

v, = (A:rki„/mo)'/2, (14) 

where mo is the mean molecular mass. With niQ « 2.3 amu for molecular clouds, 
one has w 60-y/rkin m s^' , which shows that at the typical kinetic temperature of 
10 K the nonthermal velocities are in general subsonic, and that the selected targets 
do represent the quiescent material at different distances from the core centers. 

The gas temperature in the molecular cores was estimated from the apparent line 
widths using Eq. (O. For L1498 we obtained a lower limit on the kinetic temperature 
Ticin = 7.1 ± 0.5 K (average over 8 points), which is slightly lower than T^cm = 10 K 
measured by a different method from the relative populations of the {J,K) = (2,2) 
and (1,1) levels of NH3 described by the rotational temperature T^^ BDl . For L1512 
the temperature averaged over 11 points is = 9.6 ±0.6 K, which is consistent 
with the value of 10 K. The kinetic temperatures in the L1517B core in all measured 
positions is well below 9.5 K - the value from HTI . which means that the nonthermal 
velocity dispersions of NH3 and HC3N differ significantly and that both species do 
not trace the same material. The three scanned points in the L1400K core showed 
Tkin = 8.3 ± 3.2 K - close to the expected value of 10 K (measurements of the gas 
temperature in this core were not performed in previous studies). 

We found that, in general, the spatial fluctuations of T^in do not exceed a few 
kelvin implying uniform heating and absence of the localized heat sources. The 
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Fig. 1 Example of the line-of-sight velocities of NH3 {J,K) = (1,1) (squares) and HC3N / = 
2—1 (triangles) at different radial distances along the main diagonal cut (panel a) and in the 
perpendicular direction (panel b) of the molecular core L1512. Shown are Icr error bars. For more 
detail, see (26). 

kinetic temperature tends to rise with the distance from the core center, which is in 
Une with the resuhs from BTl . 

The radial velocity profiles along the different diagonal cuts toward the selected 
targets are shown in Figs. 3 and 4 in 1261 . The diagonal cut in L1498 exhibits coher- 
ently changing velocities of V/„ (HC3N) and V/5, (NH3) except one point at the core 
edge. The velocity gradient is small, |VV/^;.| « 0.5 km s^' pc^'. This picture coin- 
cides with the previously obtained results based on observations of CO, CS, N2H+ 
and NH3 in this core and was interpreted as an inward flow 1221 |4T1 . Taken together, 
all available observations classify L1498 as one of the most quiet molecular cores. 
Thus, we can expect that the Doppler noise (irregular random shifts in the radial 
velocities between different transitions) is minimal in this core. 

In L15 12, the y/„-(NH3) and V/„-(HC3N) distributions are almost parallel (Fig.[T]i- 
The same kinematic picture was obtained for this core in l22l from observations 
of CS and NH3 lines and interpreted as a simple rotation around the center The 
velocity gradients derived from both NH3 and HC3N lines are similar, VV/^,. « 1.5 
km s^' pc^', and consistent with the gradient based on N2H+ measurements in p]. 
This means that NH3, HC3N, and N2H+ trace the same gas, so the Doppler shifts 
AV between them should be insignificant. 

In the core L1517B, which is known to be very compact ll22l . we only observed 
the central 30" x 30" where the velocities of NH3 and HC3N along two perpen- 
dicular cuts do not change much: VV),r(NH3) w 0.3 km s"' pc"', VV/.5r(HC3N) 
« 0.8 km s"' pc"^ and in the perpendicular direction Vy/ir(NH3) w VV/s,-(HC3N) 
« —0.8 km s^' pc^'. However, a wider area observation 80" x 80") of this core 
revealed an outward gas motion at the core periphery with a higher velocity gradi- 
ent, VV/^;.(N2H+) sa 1 . 1 km s^' pc^' BTl . which is consistent with earlier results on 
NH3 observations lT2l . This can lead to an additional shift in the radial velocity of 
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HC3N line since, in general, the C-bearing molecules also trace a lower density gas 
component (« < 10^ cm~^ ) in the envelope of the molecular core. A higher nonther- 
mal velocity dispersion of the HC3N line than for NH3 has already been mentioned 
above in regard to the temperature measurements in this core. 

In L1400K we only observed three positions. Both molecules trace the same gra- 
dient of Vy,„. w 1.9 km s"' pc"', which is in Hne with VV,,,(N2H+) = 1.8 ± 0.1 
km s^' pc^' derived in i4J. The mapping in different molecular lines in BTI re- 
vealed that L1400K deviates significantly from spherical symmetry and exhibits 
quite a complex kinematic structure. In particular, the distributions of N2H+ and 
NH3 do not coincide: N2H+ has an additional component to the west of the center. 
This explains why [6| reported a blue-ward skewness for this core, 9 = — 0.42± 
0.10, in the the N2H+ (1-0) hfs profiles, whereas in our observations the NH3(1,1) 
hfs transitions are fully symmetric: at the central position the skewness of the re- 
solved and single hfs component Fi,F=Q,^^l,lj of NH3 is = —0.1 ±0.3. 

Our current measurements show very similar velocity shifts (ziy)M — 25.8±1.7 
m s^' and 28. 0± 1.8 m s^' (M-estimates) for, respectively, the cores L1498 and 
L1512 where the Doppler noise reaches minimal levels. A higher shift {AV)m = 
46.9 ± 3.3 m s^' is observed in the L1517B core - again in accord with its revealed 
kinematic structure, which allows us to expect a higher radial velocity for the HC3N 
line. On the other hand, a lower value {AV)m = 8.5 ± 3.4 m s^' in L1400K may 
come from the irregular kinematic structure of the core center, which could increase 
the radial velocity of the NH3 line. 

Thus, as reference velocity offset we chose the most robust M-estimate of the 
mean value from the L1498 and L1512 cores: {AV)m ~ 26.9 ± 1.2stat m s^' . Taking 
into account that the uncertainty of the HC3N(2-1) rest frequency is about 3 m s^' , 
whereas that of NH3(1,1) is less than 1 m s^' , we finally have {AV)y^ = 26.9 ± 
1 .2stat ± 3 .Osys m s^ ' . When it is interpreted in terms of the electron-to-proton mass 
ratio variation, as defined in Eq.©, this velocity offset provides Ajj./fi — (26 ± 

Istat i 3sys) X 10 



6 Discussion and conclusions 

In two molecular cores with the lowest Doppler noise L1498 and L1512, we regis- 
ter very close values of the velocity offset AV ^21 m s^' between the rotational 
transition HC3N(2-1) and the inversion transition NH3(1,1). These values coincide 
with the most accurate estimate obtained from the Effelsberg dataset on 12 molec- 
ular clouds in the Milky Way |25J. Two other cores, L1517B and L1400K, exhibit 
velocity shifts that are either higher 47 m s^' in L1517B) or lower (~ 9 m s^' 
in L1400K) than the mean value, but the positive (L1517B) and negative (L1400K) 
deflections from the mean can be explained from the observed kinematics in these 
cores. 

Of course, simultaneous observations of the NH3(1,1) inverse transition and rota- 
tional transitions of some N-bearing molecules, such as N2H+(l-0) and N2D+(l-0) 
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would give a more accurate test. The main obstacle to this way is that the labora- 
tory frequencies of N2H+ and N2D+ are known with accuracies not better than 14 
m s^' II25I . Using the N2H+ rest frequency from the Cologne Database for Molec- 
ular Spectroscopy (CDMS) ||3T1 and observing with the Nobeyama 45-m telescope, 
we obtained a velocity shift between N2H+ andNHa of 23.0 ±3.4 m s"' in L1498, 
24.5±4.3 m s^^ in L1512, and 21.0±5.1 m s"' in L1517B |25l. We note that 
similar shifts are indicated for the central parts of L1498 and L1517B in Fig. 11 in 
II41I . where these targets were observed with the 30-m IRAM telescope also using 
an N2H+ rest frequency which is very close to the CDMS value. 

Obviously, for more definite conclusions, new laboratory measurements of the 
rest frequencies and new observations involving other targets and other molecu- 
lar transitions with different sensitivity coefficients Q are required. It has already 
been suggested to measure A -doublet lines of the light diatomic molecules OH and 
CH lITSl . microwave inversion-rotational transitions in the partly deuterated am- 
monia NH2D and ND2H IIT9I . low-laying rotational transitions in '^CO and the 
fine-structure transitions in Cl IZTl . and tunneling and rotation transitions in the 
hydronium ion H3O+ ll20l . The fourth opportunity is of particular interest since the 
rest-frame frequencies of H3O+ transitions are very sensitive to the variation of /i, 
and their sensitivity coefficients have different signs. For example, the two lowest 
frequency transitions Tr- = 1 — > 2 and 7a: = 3 J ^ '^2 °f para-HsO^ at, respec- 
tively, 307 and 364 GHz have AQ = Qt,q-i — Q^eA = 14.7, which is 4 times larger than 
AQ = 3.46 from the ammonia method. This means that the offset ZiV ~ 27 m s^' , 
detected with the ammonia method, should correspond to a relative velocity shift be- 
tween these transitions, AV ^V^^ — V^oj, of about 100 m s"' if Aji/ji ^26 x 10"^. 
We consider the hydronium method 1201 as an important independent test of the 
Afi/fi value in the Milky Way. 

To conclude, we note that in cold molecular cores with low ionization degrees 
(xe ~ 10~^ — 10^^) frequency shifts caused by external electric and magnetic fields 
and by the cosmic black body radiation-induced Stark effect are less or about 1 
m s^' and cannot affect the revealed nonzero velocity offset between the rotational 
and inversion transitions in the ammonia method. Detailed calculations of these 
effects are given in [26] . 

Our current results tentatively support the hypothesis that the fundamental physi- 
cal constant - the electron-to-proton mass ratio - may differ in low-density environ- 
ments from its terrestrial value. This may be the consequence of the chameleon-like 
scalar field. However, new laboratory measurements of the molecular rest frequen- 
cies and new observations involving other molecular transitions and other targets 
are required to reach more definite conclusions. 
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